MOND habitats within the solar system 
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Modified Newtonian Dynamics (MOND) is an interesting alternative to dark matter in extra- 
galactic systems. We here examine the possibility that mild or even strong MOND behavior may 
become evident well inside the solar system, in particular near saddle points of the total gravi- 
tational potential. Whereas in Newtonian theory tidal stresses are finite at saddle points, they 
are expected to diverge in MOND, and to remain distinctly large inside a sizeable oblate ellipsoid 
around the saddle point. We work out the MOND effects using the nonrelativistic limit of the TeVeS 
theory, both in the perturbative nearly Newtonian regime and in the deep MOND regime. While 
strong MOND behavior would be a spectacular "backyard" vindication of the theory, pinpointing 
the MOND-bubbles in the setting of the realistic solar system may be difficult. Space missions, such 
as the LISA Pathfinder, equipped with sensitive accelerometers, may be able to explore the larger 
perturbative region. 



I. INTRODUCTION 

MOND is a scheme for explaining extragalactic phe- 
nomenology without invoking dark matter. It has been 
very successful in this job despite its rather rudimentary 
form 0. In the Lagrangian formulation of MOND p\ 
the physical gravitational potential <&, which gives test 
particle acceleration by a = — V$, is determined by the 
modified Poisson equation 



V • [^(|V$|/ao)V$] 47rGp, 



(1) 



where p is the baryonic mass density, ao ~ 10" 
is Milgrom's characteristic acceleration, and the function 
ji{x) is required to approximate its argument for a; <C 1 
and to approach unity for x 3> 1. The form 



/i(a;) = x{l + xY 



(2) 



has been quite successful in modelling galaxy rotation 
curves without invoking dark matter. The theory en- 
capsulated in Eq. ^ has recently been reformulated as a 
consistent covariant gravitation theory named TeVeS H . 
Alternatives to this theory have been considered [HIEQ' 
but they shall not be employed in the present work. 

Are MOND effects of importance in the solar system 
(henceforth SS)? Milgrom was the first to consider the 
effects of MOND on the properties of long period comets 
originating in the Oort cloud 0. It was later observed 
that relativistic theories with a MOND limit can eas- 
ily predict anomalously large perihelion precessions of 
the planets!^ II]- With the discovery of the "Pioneer 
anomaly" y] much speculation was directed towards a 
possible MONDian origin of the effect ^3 ■ Relativistic 
theories with MOND phenomenology tend to produce a 
radial drift of the Kepler constant in the SS in the same 
sense as would correspond to the claimed Pioneer effect, 
though not always of the claimed magnitude S IE ^| ■ 
The persistence of Pioneer-type effects in a variety of 



scalar-tensor theories of MOND, and the hurdles faced 
by such theories from precision SS tests has been empha- 
sized by Sanders . 

In this paper we search for other sites deep inside the 
SS where strong MOND behavior might put the MOND 
phenomenon at the reach of spacecraft measurements. 
Strong MOND behavior is triggered by a low gradient 
in the total Newtonian potential $jv (the deep MOND 
regime is that where |V$| <C ao). Two apparent candi- 
dates for strong MOND regions fail this criterion. Most 
obviously we have gravitational perturbations, such as 
those accounting for the non-relativistic component of 
the perihelion of Mercury precession, or Neptune's influ- 
ence upon Uranus' orbit. Most of these have a very low 
potential gradient, and would by themselves be in the 
MOND regime. However the gradient of the total $Aris 
not small, so their effect falls in the Newtonian regime. 
Issues like the stability of the SS or its detailed dynamics 
are not expected to be appreciably different in MOND. 

The Lagrange points are another apparent possibility 
for strong MOND regions. They are the five stationary 
points of the two-body dynamics; for example LI is the 
point between the Earth and the Sun where a test mass 
would be in inertial motion, moving neither towards the 
Sun nor the Earth. Each Lagrangian point orbits the Sun 
with the same frequency as the Earth, so the gradient of 
at it must cancel the corresponding centrifugal accel- 
eration, and is thus not especially small. This does not 
mean, as we shall see, that perturbative effects around 
these points are not present; however strong MONDian 
behavior is certainly not expected. 

By contrast, the saddle (or extremum) point (hence- 
forth SP) of <I>Ar between two gravitating bodies is evi- 
dently in the deep MOND regime, since V^jv — there. 
One such point exists between any two gravitating bod- 
ies, potentially providing a testing ground for strong 
MONDian behavior. SPs are not inertial, but may be 
visited by free-falling test bodies; they are encased by 
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small "bubbles" within which strong MOND effects are 
expected. 

In what follows we study the structure of the gravi- 
tational field in increasingly smaller neighbourhoods of 
the SP for a binary mass system. Principally we look 
at the encasing Newtonian region, at the enclosed quasi- 
Newtonian sector and at the deep MOND region at the 
heart of the bubble. The strongly nonlinear character of 
the MOND phenomenon makes the analysis complicated, 
and a variety of analytical as well as numerical strategies 
have been utilized. 

The plan of this paper is as follows. In Sec. ^ we lay 
down the framework for dealing with the two (and many) 
body problem in MOND. In Sec. Illll we use heuristic ar- 
guments to find the location and extent of the principal 
MOND bubble for a binary system. One of the approxi- 
mations thus made is replaced by a more exact treatment 
in Sec. II VI In Sec.Elwe perturbatively calculate the grav- 
itational field in the quasi-Newtonian region still far from 
the SP point. In Sec. lVIl a combination of numerical and 
analytical approaches is used to deduce the gravitational 
field very near the SP where MOND is dominant. In 
Section rvTrl se discuss a number of complications to the 
above stylized treatment that arise from the many-body 
nature of the real SS. The prospect of a direct test using 
the LISA Pathfinder project [ij is briefly discussed in 
Sec. IVIIII Issues connected with the behavior of gravity 
in the spacecraft's frame are elucidated in Sec. IIXI We 
conclude in Sec. Q with a statement on how our work 
might help settling the controversy between dark matter 
and MOND. 



II. THE FRAMEWORK OF THE TWO-BODY 
PROBLEM 

Where needed we use units with c = 1. We base 
the analysis on the non-relativistic limit of TeVeS 0; 
AQUAL and other Lagrangian formulations of MOND 
have similar form and many of our conclusions may ap- 
ply to aU of them. In TeVeS the MOND behaviour is 
driven by a dynamical (and dimensionless) scalar field 4> 
such that the physical potential $ in which a body falls is 
given by $ = $Ar -|- 0, where $Ar is the usual Newtonian 
potential (inferred from the metric component goo)- In 
the non-relativistic regime is governed by the equation 



V • [Ai(fcZ2(V0)2)V0] = kGp 



(3) 



where fc is a coupling constant and / is a length scale 
which determines the Milgrom acceleration by 



3fc 



oo 



47r/ 



lO-^V s-2 



(4) 



(we are setting S, as defined in to unity; thus we ig- 
nore the slight renormalization of the gravitational con- 
stant in TeVeS so that here Gn = G). In Eq. © /i 
is a free function not to be confused with Milgrom's jl. 



Ref. 1^ proposed a particular form for it. The deep 
MOND regime is signalled by the low gradient of the 
scalar field 0; in this regime 



fc |V0| 
At: flo 



(5) 



For strong gradients the /i proposed in Ref. 4] grows 
crudely as (| V0|/ao)^/'^. This has the effect of suppress- 
ing the contribution of V0 to V<I> thus bringing in the 
Newtonian regime. In spherically symmetric systems 
TeVeS with any /i satisfying Eq. Q goes over into the 
Lagrangian MOND theory with Milgrom's /i given by 
/2 = (1 + fc/47r/^)^^. Although this point is not well ex- 
plored, it is quite possible that in less symmetric systems 
TeVeS does not go over to an exactly MOND behavior. 
For this reason we base this paper on the non-relativistic 
limit of TeVeS, and not on Lagrangian MOND. 

We need to solve Eq. Q for a two-body source, but 
that equation is non-linear so the fields due to each 
body do not superpose. However any non-linear equation 
may be formally linearized by an appropriate change of 
variables. Here this is 



V0 



(6) 



(see [13, where this technique was first suggested). We 
may then add the u due to each source (which is the 
Newtonian acceleration) and invert the total u at a given 
point to find V0. It is essential that the sum of all sources 
be performed before inverting to find V0. 

This algorithm may be applied to any number of com- 
ponents. But note that even if a term in the sum is in 
the MOND regime, the overall system is not, unless the 
total |u| is much smaller than oq. (It is because of this 
feature that the gravitational perturbations in the SS are 
non-MONDian.) However, it is also possible to have two 
components with fields not in the MOND regime such 
that their common field is MONDian in some region. Ex- 
amples are the SPs in the gravitational potential of two 
bodies to be studied in this paper. 

The only complication with the above technique is that 
u is generally not curl-free; indeed it is rather the vector 
u/^ which is curl-free. Thus the full set of equations for 
u is 



V • u = -AirGp 
V A - = 0. 



(7) 
(8) 



The first equation tells us that u equals the Newtonian 
acceleration F^^' = — V<&Ar up to a curl, that is, there 
must exist a vector field h such that 



u = f(^' + V a h. 



(9) 



The second equation fixes the h (up to a gradient). This 
operation can only be performed upon the total u, once 
again stressing the intrinsic non-linearity of the theory. It 
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can be shown that the curl term vanishes in a spherically 
symmetric situation, or in the quasi-Newtonian regime 
far away from the source 0, Q . Near the SPs neither 
of these conditions is satisfied and we have to evaluate 
V A h. However, before plunging into the full problem, 
let us provide some orientation 



III. HEURISTICS OF THE MOND BUBBLES 

Consider two bodies at distance R with masses M and 
TO, with M ^ TO, so that the system's center of mass 
may be taken to coincide with the heavier body. To be 
definite we call them the Sun and the Earth, but we shall 
explore other couples later. Along the line linking them 
(the z axis), the Newtonian acceleration is 



GM 



Gr 



(10) 



where r is the distance from the Sun and is the unit 
vector in the direction Sun to Earth. The SP of the 
Newtonian potential resides where F^^^ = 0, i.e. at 



r = rs ^ R I - 




(11) 



Around this point F*^^) increases linearly as it passes 
through zero, that is 



fW « A(f-r,)e, 



where 



A 




(12) 



(13) 



is the tidal stress at the SP along the Sun-Earth direction. 
The full tidal stress matrix is easy to compute. Let us use 
cylindrical coordinates centered at the SP, with the z-axis 
pointing along the Sun-Earth direction, so that we have 
dF^^^/dz = A and dP'^^'^ /dz = 0. From the further 
condition that the divergence must be zero (outside the 
Sun and Earth) we have 



fW =A(ze,--ge,). (14) 

The region around the SP is obviously in the deep 
MOND regime (|F(^)| < ao). Thus regardless of the 
model adopted for /z, we have just on the basis of Eq. © 
that 



u ; 



ao 



-V0 



V Ah. 



(15) 



We may estimate the physical acceleration acting on 
an object by ignoring the curl term here. We see that 
|V(/)| < |F(^)|, so that 



F = -V$ 



2 



1/4 



(16) 



In contrast to the Newtonian theory, the tidal stresses 
here diverge at the SP. This may be understood by ap- 
plying the rule of thumb that in the deep MOND regime 
the square root of the Newtonian acceleration gives the 
physical acceleration. According to Eq. (|12() Newtonian 
acceleration increases linearly along the line Sun-Earth, 
so the physical acceleration in the deep MOND regime is 
of the form ± ^ao|r — Ts |, which has infinite derivative 
at Vs. 

The tidal stresses are expected to remain anomalously 
high, and gravity to remain in the deep MOND regime 
(with |V0| ^ |V$7v|), within a very small oblate ellip- 
soidal region around the SP. The size of it can be esti- 
mated from the condition |F*^^^| = ao which translates 
into a major semi-axis (in the g direction) of size 



6f 



(17) 



where a„j is the acceleration of the smaller mass m. The 
semi- minor axis (aligned with the z axis), has half this 
size. For the Sun-Earth system this is very small, Sf « 
4.4 m, but it gets larger for the outer planets, for example 
for the Sun- Jupiter system we have i5f « 11 Km. The 
Earth-Moon system gives i5f « 1.6 m, but this is actually 
the best ratio 6f/R of the three examples. 

There is a much larger intermediate region where there 
are significant perturbative corrections to Newtonian the- 
ory, but where deep MOND behavior is not yet in evi- 
dence. However this transition region is very model de- 
pendent. For the model introduced in Q Milgrom's fl 
can be estimated in the quasi-Newtonian region through 
formula (69) there (which formula, however, is rigorous 
only in the spherically symmetric case): 



A* 



1 



k3 p2 



(18) 



Let us take F « i^(^) and use Eq. (HH). We see that 
10~^ departures from Newtonian gravity occur within a 
semi-major axis of size 



8007r3/2 
fc3/2 ~A 



(19) 



Using k « 0.03 (as suggested in 4]) this is Af 



1,900 Km, Af « 4.7 x 10^ Km, and Af 



700 Km, 



respectively, for the three examples described above. 

Of course the calculation in this section is very ide- 
alized. Orbits are not circular, the center of mass does 
not coincide with the Sun center, and rather than two 
objects we have a multitude of competing influences. In 
Section rVIII we consider more realistic situations, examin- 
ing also the surprising effect of the extra solar component 
(negligible in Newtonian theory but not in MOND). Be- 
fore that, however, we must reconsider the effect of the 
neglected curl term. We do this in Sections [l VI VII While 
the qualitative aspects of the present section survive the 
introduction of the curl term, we shall find that the quan- 
titative conclusions are significantly modified. 
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IV. ACCOUNTING FOR THE CURL TERM 

By carrying out the curl in Eq. ^ we get 

Vln/iAu-VAu = 0, (20) 

while squaring Eq. JHjl gives 

= (4Vfc)V|V</)p. (21) 

In TeVeS /i = /i(fc|V(/)|/ao); thus k^u^ /ao^ is a function 
of /i only. Defining the dimensionless quantity 

K = ^\TLU^/^\TL^l, (22) 

we get Vln/i = k^^Vu'^/m^ so that Eq. H20() becomes 

KU^VAu + uAVu^ = 0. (23) 

In systems with spherical, cylindrical or planar sym- 
metry, u is necessarily coUinear with V|up. Then V A u 
must vanish everywhere (since n would be expected to 
vanish only at isolated points). This agrees with the 
findings of Refs. 0-0] that /zV(/) and /xV$ are both curl 
free in such situations. When the spatial symmetry is 
lower or nonexistent, the second term in Eq. (|23|l will 
not generally vanish, and will be of order |up/L where 
L denotes the scale on which quantities vary. Thus if in 
a region k ^ 1, we would expect |V A u| to be much 
smaller than its expected magnitude |u|/L; this signals 
the quasi-Newtonian regime where u is nearly curl-free. 

In TeVeS the manner of transition between the deep 
MOND and Newtonian regimes is dependent upon the 
form of /i. The form proposed in Ref. 0| is quite difficult 
to work with in our context. We shall thus replace it by 
the implicit expression 



P(x) 



k |V(/)| 
47r ao 



(24) 



which satisfies the limit (0). A simple calculation (see 
Ref. 0) then shows that Milgrom's //, defined by Eq. 
is here given parametrically by 



X = C(l + C) [1 - (fcC/47r) 



41-1/2 



(25) 
(26) 



which satisfies the MOND requirements that pi ^ x for 
a; <C 1 and fl —>■ 1 as x ^ oo. Fig.Hcompares our jl with 
the "simple" of Refs. Q and Q 

Eliminating |V(/)|/ao between Eqs. H21|l and 124|l gives 



2567r4 



Differentiating the logarithm of this we calculate that 



(27) 



647r4 ao^ 



(28) 




FIG. 1: Milgrom's jl function (upper curve) as used in many 
extragalactic applications (Eq. J^J) and as defined in this pa- 
per with k = 0.03 by Eqs. 



In terms of the dimensionless vector field 

fc2 u 



Uee 



we may thus cast Eqs. H23|l and ((Tj) into the form 



(29) 



V • U = (30) 
4(l + t/2)C/2VAU + UA = 0. (31) 

where we have dropped the source of the first since we 
are interested only in the region near the SP. This pair 
of exact equations for one dimensionless vector is central 
to our study. 

Once U is solved for we can recover V0 by combining 
Eqs. ®, (123 and HHj): 



V0 : 



^™°(1 + C/^)V4-^ 



k 



C/1/2 



(32) 



As remarked earlier, the condition k ^ 1 brings in the 
Newtonian limit. Now k 3> 1 is equivalent toll ^ 1. Ob- 
viously in this case — V0 w (47rao/A:)U = (A:/47r)u which 
tells us by Eq. ^ that /i ~ 1, indeed the Newtonian limit 
(the same is obvious from Eq. (HSJ). 



V. THE QUASI-NEWTONIAN REGION 

At this point we go over to spherical polar coordinates 
(r, '0, 0) with origin at the SP; accordingly 



z — r cos "0; g ~ r sin ip. 
So, for example, Eq. H14I) takes the form 
F(^) = ArN, 



where 



N(V') 

Nr 



-[l + 3cos(2V')] 
--sin(2^/>). 



(33) 

(34) 

(35) 
(36) 

(37) 
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FIG. 2: The flow of Uo around the SP (at the origin) in a 
plane containing the symmetry (z) axis; for clarity all vectors 
have been linearly rescaled. 



We define the quasi-Newtonian region as that where 
is of order 1 or larger so that the factor 1 + cannot be 
ignored in Eq. The region's size may be estimated 

by dropping the curl term in Q (an approximation to be 
justified a posteriori) and finding the solution to = 1 
using l|14|) and 129|) (in the Newtonian region u = F^^^). 
This leads to the ellipsoid: 



■(cos2^+-sin2^j =r2^ (^5^ 



(38) 



Eq. H31|l tells us that well outside of this ellipsoid the 
curl is suppressed by a factor of with respect to 
F^-'^). As we show below U is then neatly separated 
into a Newtonian component Uq (carrying the divergence 
predicted by iQ and depicted in Fig.|2Jl and a "magnetic" 
component U2. By definition U2 is solenoidal and to 
leading order is sourced purely by Uq. Specifically the 
dynamics is approximated by 



U 

Uo 

V-U2 

V AU2 



Uo 

r 







U2 

N(V) 



Uo AVlUol 



With the notation 

U2 = Urer 

Eqs. (|^ and become 



j,2 Qj. 



1 



d 



r sin ^ d^p 



4|Uo| 



(sin-f/iC/^) 



(39) 
(40) 

(41) 
(42) 

(43) 

0, (44) 
(45) 




FIG. 3: The angular profile functions F (solid) and G 
(dashed) giving the direction of the "magnetic" field B in 
the quasi-Newtonian region; for clarity G has been multiplied 
by 3. 



with 



cos ip sin ip 



9 / I sin^ ijj 
COS'^ ip -\ 3-^ 



12 sin 2^- 
(5 + 3cos2V')2 



(46) 



The form of Eqs. and suggests that both Ur 
and behave as 1/r. Accordingly we recast Eq. H43|) 
as the ansatz 

U2 = ^B(^) = ^ (^^(V^)e, + G(V')e^) , (47) 

where the r dependence has been fully factored out. With 
this ansatz Eq. H45|l collapses into 



F' ^-s^ 



12 sin 2^- 
(5 + 3cos27/')2 ' 



with solution 



F 



A, 



5-1-3 cos 2ip 
where A is a constant. Eq. H44() now becomes 
1 d 



sin %p dip 



which integrates to 



G sin -0 



(sinV'G) = 0, 



F sin Ip d-p + B 



(48) 



(49) 



(50) 



(51) 



where B is another constant. Performing the integral 
gives 



Gsimp 



ta.n-^{V3 - 2 tan |) + ta,n-^{VS + 2 tan |) 

+Acosip + B. (52) 

To determine A and B we must discuss boundary con- 
ditions. According to Milgrom for the system H30() - 
(|31|l the normal component of u (or U) must vanish on all 
boundaries. Parts of the symmetry axis (■0 = as well as 
■0 = tt) are evidently a boundary of the quasi-Newtonian 
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FIG. 4: The flow of U2 in a plane containing the z axis; 
coordinates are in units of ro . For clarity the solution was cut 
off at r = ro (so as to avoid a divergence at the origin). 



region; it is obvious that vanishes on both North 
and South parts of it, where it is the normal component. 
Thus since Uq satisfies the boundary condition on the 
relevant pieces of the axis, so must U2 . Accordingly we 
must require G(V' = 0) = G'('(/' = tt) = 0, from which 
follows that 



3\/3 ■ 



(53) 



The solutions F and G are plotted in Fig. |31 We find 
that G(7r/2) = as well. Thus on the symmetry plane 
(-0 = 7r/2) U is coUinear with the axis. 

What about the rest of the boundary? We see from 
Eq. (|47|l that U2 ^ as r — > 00. Thus at large r our U 
merges with Uo which we know to be the limiting form 
of the Newtonian field as we approach the SP. It follows 
that our solution automatically fulfills the boundary con- 
ditions at large r. The inward part of the boundary of the 
quasi-Newtonian region adjoins the intermediate MOND 
region, where MOND effects are no longer small. Fortu- 
nately there is no need for us to set boundary conditions 
there; rather, the solution just described serves to set 
boundary conditions for the intermediate MOND region. 

We conclude that a SP far away from the strong 
MOND bubble is characterized by a Newtonian compo- 
nent proportional to r together with a magnetic-like per- 
turbation that falls off like 1/r. The full physical effects 
in this regime may be appreciated by combining (|32|l with 
(I40|l and (|47|) . We find that the extra acceleration felt by 
test particles is 



Up 

ml 



(54) 



The first contribution, call it JFq is of fully Newtonian 
form, and just serves to renormalize the gravitational 
constant, as discussed in Ijl. The second term was also 



derived in (c.f. Eq. (69) of ^\) and is 



(5Fi ^ 



fc3 p{N)2 



k r 5 -I- 3 cos(2'(/') 



(55) 



What we have just shown is that to these two terms one 
should add the magnetic-like contribution 



^F. = i!;^^r£B(^), 



(56) 



which is of the same order of magnitude as S¥\. Apart 
from the prefactor -iTTao/k, this term is just what was 
plotted in Fig. 0] In Fig. [S] we plot the angular profile 
B(-0) + 2[5 + 3 008(2-0)] "^N(V') of the total correction to 
the acceleration after renormalization of G. The plotted 
field is to be divided by r (and multiplied by 47raoro/fc) 
to obtain the extra acceleration felt by test particles in 
the quasi-Newtonian region. 

How do these results affect the naive expectations of 
Section Init ' We have just shown that a full quantitative 
analysis can never neglect the "magnetic" field derived 
in this Section. In addition the border between full and 
linear MONDian behavior is determined by the condi- 
tion = 1, equivalent to ellipsoid H38() . As long as we 
stay well outside this ellipsoid we obtain results consis- 
tent with H18|l and (|19|l : however the order of magnitude 
of linear corrections outside this ellipsoid may be written 



dF 



We learn that the highest fractional correction in this 
regime is achieved close to ellipsoid (I38f) and is of order 
fc/(47r), around a 0.0025 for k w 0.03; it then falls off as 
l/r^ as we move away from the SP. Therefore, as long 
as we do not use (|18|) for fractional corrections larger 
than A:/(47r) we obtain qualitatively correct results (the 
example given in Section [llll satisfies this condition). 

The bottom line for our predictions is that the ellipsoid 
(|38|l represents both the region where the largest linear 
corrections are felt and the border for the onset of full 
MOND behavior. For the three examples considered in 
Section HTll we have 



ro 



383 Km 



Earth — Sun 



(58) 



ro « 9.65 X 10^ Km Jupiter - Sun (59) 
ro « 140 Km Earth - Moon (60) 

corresponding to ellipsoids with major semi-axis of 766 
Km (Sun-Earth), 1.93x10^ Km (Sun- Jupiter) or 280 Km 
(Earth-Moon). These are the relevant dimensions of the 
MOND bubbles. 



VI. THE DEEP MOND REGION 

By Eq. ^ the deep MOND regime (/x <C 1) entails 
w 4 or f7 < 1. Thus in Eq. (|2U we replace 1 -I- ^ 1. 
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FIG. 5: The flow of B{iP) + 2[5 + 3 cos(2V')]"^N(i/>). When 
divided by r and muhiplied by 47raoro/fc this fleld gives the 
physical acceleration beyond the Newtonian one felt by test 
particles in the quasi-Newtonian region. 



Then together with Eq. H30() this has a double symmetry, 
already noticed by Milgrom 15j. They are both invariant 
under U const, x U (rescaling), and under x Ax 
(dilation of the coordinates). The first symmetry implies 
that the normalization of U is arbitrary (of course the 
normalization is eventually fixed by taking cognizance of 
the sources of Eq. (|30l) ). The second means that a solu- 
tion whose linear scale is expanded remains a solution. 

In spherical polar coordinates Eqs. (|30|I - H31II take the 
form 



1 d 



J.2 Qj. 

r\ dr dip 



1 



d 



r sin ij} dtp 
Ur_d_ 

r dip 



(sin-0 
d 



dr/ 



0(61) 
0(62) 



For a solution of these to turn into a second solution 
upon dilatation of the coordinates {r — > A r) , it is nec- 
essary for the r dependence of both Ur and to be a 
single power. Thus we make the ansatz 



U = C 



(63) 



with C and a dimensionless constants. The power a — 2 
was chosen for notational convenience in what follows. 
Substituting in Eq. H61|l we obtain 



G" -I- ctan(i/;) G + aF 
while substitution in Eq. 1)62(1 gives 







(64) 



F- 



dtp 



2[aG-2F']{F^ + G^) =0. (65) 



These last constitute a coupled system of first order or- 
dinary differential equations for F{tp) and G{tp). 

These equations have several symmetries. F and G 
may be rescaled, that is multiplied by a constant (this 




FIG. 6: The numerically determined angular profile functions 
F and G in the deep MOND region (solid) compared with the 
Newtonian profile functions Nr and A'^^ (dotted), respectively. 



is nothing but the scale-invariance of the deep MOND 
regime). We also have the symmetry: a — *■ —a, F F, 
G — *■ —G. Finally the equations are parity-invariant: 



"0 —>■ TT — "0, 

F -> ±F, 
G ^ TG. 



(66) 



Of course this by itself does not compel the solutions 
themselves to have definite parity, i.e. F{ip) = ±_F(7r — 
ip) and G{ip) = TG(7r — tp)- However, numerically we 
find that the only regular solutions are indeed those with 
definite parity, and that these only exist for a discrete 
sequence of as: {±q;i, ±q;2, • ■ ■ }. Specifically we find 
ai = 2 and the approximate values a2 ~ 3.528, as ~ 
5.039, a4 ~ 6.545, etc. 

Seen in another way, the boundary conditions at -ip — 
(see below) justify representing as a Fourier series in 
cos^rmp) and G as a Fourier series in sin(m7/)). It is only 
for the mentioned special as that even and odd m modes 
decouple, so that we can have a solution that is a series 
in only odd or only even m. For other values of a the 
solutions mix even and odd m, but are singular at ip ^ n. 

We further found numerically that each of the regular 
solutions is dominated by a single Fourier mode, i.e. it is 
of the form F f^i Fq + Fn cos mp and G « G„ sin mp for a 
given n. The solution corresponding to n = 1 is obtained 
for |ai| = 2 and is exactly F — acostp and G — ^asintp 
(a is a constant) for a = ±2, respectively. We have been 
unable to find analytic expressions for other coefficients 
Fn and G„, but have determined them numerically. We 
now select the relevant solution by imposing appropriate 
boundary conditions for our problem. 

As in Sec. the boundary condition that the normal 
component of U vanish requires that we take Giip = 0) = 
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G{ip = tt) = 0. Because C can still be adjusted, we loose 
no generality in requiring the corresponding boundary 
condition F{ijj = 0) = F{ip = n) — 1. For were we to 
demand F{'>jj = 0) ^ F{ip = tt), we would thereby intro- 
duce a jump in U across the plane = 'k/2 for which 
there is no physical reason. Our choice of boundary con- 
ditions immediately selects a solution with definite par- 
ity, which as mentioned earlier, are the only nonsingular 
ones. Regarding boundary conditions at large r, we know 
that there must be a match with the field in the quasi- 
Newtonian region. This naturally selects the particular 
solution with n — 2, since the quasi- Newtonian solution 
Uo has components with angular profiles of form cos 2ip 
or sin 2-0. This logic still does not prefer positive over 
negative a. But to avoid a singularity at the origin (see 
Eq. (|(j3|l 'l we should select the solution with positive a, 
namely that for a ~ 3.528. 

The functions F and G obtained for this a are plotted 
in Fig. El These graphs are approximated at the level of 
1 % by the formulae 

F(V') = 0.2442 + 0.7246 cos(2V') + 0.0472 cos(40), 
G(0) = -0.8334 sin(2V') - 0.0368 sin(40). (67) 

For comparison Fig. El plots also Nr and of 
Eqs. (|35|I -H37 |I . We see that the angular profile of the 
deep MOND U (whose flow is plotted in Fig. [7|) is quite 
similar to that of the Newtonian Uq (Eq. and Fig.jSJ. 
Of course, the radial dependences of the two are quite dif- 
ferent. Now as mentioned earlier, in the absence of any 
mention of the sources in Eqs. H30|) - (|31|) . it is not possi- 
ble to determine the normalization of U. However, we 
may estimate C in Eq. H63|l as follows. Given the simi- 
larity of the angular profiles wc may suppose that were 
we to extend the deep MOND U of Eq. ^ to the inner 
boundary of the Newtonian region at r = ro, we should 
obtain Uq. This requires that C = 1 and we adopt this 
value. 

We conclude that taking the curl term into account in 
the deep MOND regime once again vindicates qualita- 
tively the simplified arguments of Section IlIII but intro- 
duces substantial quantitative novelties. Using we 
find that the extra physical force is now 



Anao U 
k C/i/2 



(68) 



If we define D as the angular profile in the deep MOND 
regime (which as we have seen is very close to N) then 



SF 



47r an 



D 



(69) 



For a < 4 (a condition satisfied by our solution) , the tidal 
stresses associated with this field, i.e its spatial deriva- 
tives, diverge at the saddle, as predicted in Section ITTll 
However the divergence is softer than in Eq. H16|) where 
the curl term was ignored (that solution corresponds to 
a = 3, which is an unphysical value as we have seen). 
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FIG. 7: The flow of the fleld U in the deep MOND regime 
plotted with linear scale in units of ro and assuming C — 1. 



Rearranging (|69|l with the aid of the definition of rg 
and Eq.l^we find that the continuation of formula l|57|) 
in the deep MOND regime is 



SF 



47r V r 



47r V r 



-0.24 



(70) 



Hence the fractional correction to Newtonian gravity, 
which equals fc/(47r) « 0.0025 at the ellipsoid con- 
tinues to grow in the strong MOND regime as we ap- 
proach the saddle. (Were we to ignore the curl term, 
in which case SF/F^^'^ ^ this growth would be 

steeper.) One implication of the growth is that the 
(f) force overtakes i^(^) in a much smaller inner region 
than naively expected (cf. formulas p7|l V Specifically 
SF w F(^) at 



A\Att 



(71) 



This is smaller than (|17|l by a factor of 10"^, and is es- 
sentially microscopic except for the Jupiter-Sun system. 
The value of F^^^ when it becomes subdominant is not 
flo as naively expected; it is also smaller by a factor of 
10-6. 

In summary, the full analysis reveals that there is a 
very large region (given by the ellipsoid (|38|l ) inside which 
fuU MONDian effects are present. The fractional MON- 
Dian corrections to gravity in this region exceed fc/(47r) 
and are therefore significant. However the MOND field 
only dominates the Newtonian field, i.e. the fractional 
correction becomes larger than unity, in a region far too 
small to be observable. 
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VII. THE REALISTIC SOLAR SYSTEM 

The results of Sees . IIIII IVl and IVll can be used to show 
that the SP location for a pair of masses as determined by 
pure Newtonian gravity (F^^^ = — V<f>Ar = 0) coincides 
with that determined by full TeVeS (-V($Ar + <f)) = 0). 
In the calculations in Sec. lVII the origin r = is the point 
where V0 = 0, and the field configuration of W(f>, or its 
surrogate U, in North and South hemispheres are reflec- 
tions of each other (see Fig. [7^ . This configuration acts 
as a boundary condition for U in the quasi- Newtonian re- 
gion (treated in Sec. 01. Accordingly, we expect not only 
the "magnetic" part U2, but also the Uq, which serves 
as background for U2's equation H42|l . to reflect the men- 
tioned symmetry, and for the null points of both these 
fields to coincide with that of U of the deep MOND re- 
gion. Now as we move outward from the quasi-Newtonian 
region, U becomes dominated by Uq which is the pure 
Newtonian field. Hence the SP determined by that field 
(see Sec. IIII|) coincides with that determined by the full 
MOND field. 

The results presented so far are "model calculations" , 
valid under a number of simplifying assumptions not sat- 
isfied by the real SS. For example, orbits are elliptic, not 
circular; the barycenter of the system does not coincide 
with the center of M; we have a many-body problem, 
not a two-body problem; etc. To leading order these 
complications do not change the anomalous effects pre- 
dicted by MOND around SPs or the size of the regions 
where they are felt. They do complicate the issue of lo- 
cating "MOND bubbles" , but since their centers coincide 
with the SPs of the Newtonian potential, this is in fact a 
Newtonian physics problem, independent of MOND dy- 
namics. 

For example the SS barycenter is dominated by the 
Sun-Jupiter pair and lives just outside the solar surface, 
rotating with a period of approximately 11 years. But 
even this is a crude approximation: the relative position 
of the Sun and any planet depends on the configuration 
of the entire SS, and is chaotic. The same may be said for 
the location of the SP between the Sun and that planet. 
However, with empirical inputs and a numerical Newto- 
nian code we can determine the location of the full set 
of SPs, and even predict where they will be within a few 
years [Ig • Not only are these details in the realm of New- 
tonian physics, but they do not affect our conclusions on 
MOND effects around SPs, as long as we anchor our so- 
lutions to wherever Newtonian theory predicts the SPs 
to be. Indeed we only need the result H14I) as a boundary 
condition for our MONDian calculations. 



effect of the Moon, the Earth-Sun SP is predicted to be 
well inside the Moon's orbit, so we cannot decouple the 
two systems. This induces dramatic qualitative changes 
in the location of both SPs. 

Let us assume that within the Moon-Earth orbit the 
Sun's field is uniform (with strength a^) and parallel to 
the Lunar orbital plane (this may be refined, but does not 
alter our point). If we consider a frame attached to the 
unperturbed Moon-Earth SP, with z pointing away from 
the Earth and x on the orbital plane, then the Newtonian 
acceleration for points on this plane is 

Fn — A(^zez xBx) aE[cos{Lut)ez — siii{ujt)ex], (72) 

where A is given by Eq. H13|) with m the Moon's mass 
and M the Earth's, uj = 27r/T, T the Moon's synodic 
period, and t ~ the time when the Sun, Earth and 
Moon are aligned (in that sequence). Thus, within these 
approximations, the actual SP (defined by Fat = 0) de- 
scribes a monthly ellipse centered on the unperturbed SP 
location, with equation 

X = 2^sin(cjt), (73) 
z = ^cos{Lut). (74) 

Thus the semimajor axis is 2aE/A^ 6.0 x 10* Km while 
the semiminor one is qe/A « 3.0 x 10* Km; these axes 
perform an annual rotation so that the smaller axis stays 
aligned with the Sun. Since the rough prediction (ignor- 
ing the Sun) places the SP about 4.3 x 10* Km away from 
the Moon, we cannot ignore the details arising from the 
three-body problem. When all effects are taken into ac- 
count, the SP may encounter the Moon's surface in its 
motion. 

Other perturbations superpose further ellipses upon 
this motion, each aligned with the direction of its source. 
In general the motion of the SP is a series of elliptic epicy- 
cles of this form due to all possible perturbations on the 
main two-body system. 

It should be pointed out that the fact that the Earth- 
Moon system falls freely in the field of the Sun does 
not alter the above arguments (see Sec. IIXII . In spite 
of the weak equivalence principle, the criterion for strong 
MOND effects is that the field F^^) calculated in the 
global frame becomes comparable to ao so that V(/) can 
make a significant and identifiable contribution to the 
overall gravitational field — V(<f>Ar + (p). Hence the posi- 
tion of the SP is to be determined, to first order, as done 
above. 



A. An example: Epicycles and the many-body 
problem 

However, these practical details do affect the planning 
of experiments, as we now illustrate. Consider the loca- 
tion of the Earth-Sun and Moon-Earth SPs. Ignoring the 



B. MOND bubbles as accelerometers 

There is a fine detail of the SP system that is purely 
MONDian: sensitivity to the extra-solar potential, or 
more precisely, to the peculiar acceleration, a^, of the SS 
barycenter. As explained before, MONDian behavior can 
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only be identified from the total potential and this must 
include the extra-SS component. In Newtonian theory 
the effect of ap passes unnoticed because nothing dra- 
matic distinguishes the SP; by contrast in MOND the 
SP is signaled by diverging tidal forces. The fact that 
we are free-falling in this field is irrelevant as just men- 
tioned. We note that ap does not shift the location of the 
Lagrange points because their definition involves balanc- 
ing inertial forces and gravitational fields. Since we are 
free-falling in the extra-galactic field, aside from making 
a tidal distortion, this acceleration does not affect the 
Lagrange points. 

The effect of Up is to superpose a further elliptical mo- 
tion onto the larger epicycles due to intra-solar perturba- 
tions. If Hp = {upz, ttpx, Upy} in the system of axes used 
above, then a similar calculation leads to 

2 

z = —{apzCoa{uj't) — apxSm{uj't)), (75) 
2 

X — —{ap2sm{uj't)+apxCos{uj't)), (76) 

y = ^(^py^ (77) 

where lu' is 2t: divided by the Moon's sidereal period. We 
see that the SP is raised off the orbital plane by apy/A, 
and describes an ellipse oriented with ap on this plane. 

How can we estimate ap? This is the "acceleration 
counterpart" to the CMB dipole (which measures the pe- 
culiar velocity Vp with respect to the cosmological frame). 
There is no simple way to estimate ap other than iden- 
tifying all components making up Vp and inferring the 
acceleration on a case by case basis. Part of Vp is due 
to our motion around the Milky Way at 217 Km/s and a 
radius r of about 26,000 ly. From these figures we infer 

a = w 1.9 X 10"^"m s"^ (78) 

which is comparable to ap. In addition there are "non- 
linear" peculiar velocity components, such as the move- 
ment of the Milky Way about the center of the local group 
and the motion of the local group toward the great at- 
tractor. These are of the order 100-200 Km s~^, and 
deriving their associated acceleration is complex. 

When all these non-linear components are added, the 
total points in the opposite direction to the CMB dipole 
(which implies a speed of roughly 300 Km s~^), so we may 
conclude that the peculiar velocity is about 600 Km s^^ 
roughly in the direction of the CMB dipole. In cosmolog- 
ical linear perturbation theory there is a simple relation 
between peculiar velocity and acceleration, namely 

ap = -QP„^^HoVp, (79) 

where Hq is the Hubble constant and f2„j is the ratio 
between the matter density p,„ and the critical density 
Pc ■ For currently popular values of these parameters 
we have Op « 1.3 x 10~^^m s~^. (Naively this places ve- 
locity perturbations in the MOND regime, but the crite- 
rion for MONDian behavior on cosmological scales should 



be derived from a MONDian counterpart of the above 
mentioned perturbation theory; this is just now becom- 
ing possible 0|.) 

We may guess that (|79|l provides a good order of mag- 
nitude estimate for the acceleration due to both linear 
and non-linear large-scale perturbations. This suggests 
that the SS peculiar acceleration is dominated by its 
motion around the center of the galaxy. However, one 
cannot discount the possibility that the SS has a signif- 
icant mass in its neighborhood, e.g. undiscovered mas- 
sive planets or even a stellar companion. This would con- 
tribute to Op. In any case, this effect has been constrained 
using timing data on accurate astronomical clocks |l9j| . 
leading to the bound Op < 5 x lO^^^m s^^. 

The shift of the SP due to ap is small, on the order of 
meters for the Earth-Sun system. But should MONDian 
behavior as predicted in this paper be discovered, the 
motion described by the SPs provides our best chance for 
a direct measurement of the peculiar acceleration; MOND 
bubbles would then function as sensitive accelerometers. 



VIII. TARGETS FOR LISA PATHFINDER 

As stated in the introduction, the MOND effects near 
the Lagrange points are expected to be weak; however 
this does not mean that they are beyond the reach of 
very sensitive equipment, such as that on board of the 
LISA Pathfinder (LPF) mission ■ Furthermore, while 
in transit to LI, the satellite may pass close enough to 
the SP to probe the quasi-Newtonian region examined 
in Section (the extreme MONDian region described in 
Section IVTl probablv requires a dedicated mission). In 
the LPF mission two proof masses are suitably shielded 
from radiation pressure and other annoyances that pre- 
vent testing gravitational physics to accuracy in the 
inner parts of the Solar system. Naturally the satellite 
itself has to bear radiation pressure, but its orbit is cor- 
rected by tracking the free falling proof masses contained 
in its inside. The sensit ivity to tidal stresses has been 
quoted as IQ-^^s-^ (see HI). 

According to Eq. H13|) . tidal stresses at the Sun-Earth 
SP are of the order A k, 4.57 x 10~^^s~^, four orders of 
magnitude larger than LPF's sensitivity. The fractional 
corrections to Newtonian gravity contained in Eqs. 155|) 
and (|56|l . and plotted in Fig. |31 have a rough order of 
magnitude given by (|57|l . The tidal stress corresponding 
to 5F is thus of order 10~^^(ro/r)^ s~^ for the illustrative 
value k = 0.03 used in this paper. Therefore LPF would 
be sensitive to these MONDian corrections if it got to 
within lOro ~ 3830 Km of the saddle. This is not overly 
demanding; the region is the size of a planet. The MOND 
effects may be even apparent while LPF is in transit to 
LI. 

In contrast, the MONDian tidal stresses felt near LI 
are far too small to be at reach of this mission. If 
denotes Ll's distance from the Sun, LI lies at i? — « 
1.5 X 10^ Km from Earth; the saddle of the Sun-Earth 
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potential is at i? — « 2.6 x 10^ Km from Earth (see 
Eq. Therefore LI is Ar = -r^ « 1.24 x 10^ Km 

away from the saddle, implying suppression of corrections 
to Newtonian gravity by a factor ^(^) ~ 2.4 x 10^^". 
By way of contrast, the Newtonian tidal stresses at LI 
are, say for the radial component, 



dF, 



(N) 



dr 



8^1 



3.17 X 10 



(80) 



with uje being the angular frequency of the Earth's orbit. 
This is only 2 orders of magnitude above experimental 
sensitivity, and so the MONDian corrections to stresses 
in the vicinity of LI are 8 orders of magnitude too small 
for the quoted instrumental sensitivity. 

However, "indirect" effects may possibly be detectable 
by LPF: effects not on its accelerometers but on its path 
(this comment may apply to other LI missions). Indeed 
MOND introduces a small shift to the location of LI and 
its surrounding orbits. Combining Eqs. H55() . (|56|l and 
(|49|l . we obtain an extra acceleration at LI with radial 
component of signed magnitude 
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An Tq 
k Ar 
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2 3V3. 



ao « -1.3x lO^^Vs"^ . 



Hence this extra acceleration predicted by MOND points 
towards the SP, i.e. away from the Sun and toward the 
Earth. In the usual calculation, the centrifugal acceler- 
ation at LI is exactly balanced by the gravitational one 
F(^). This last has absolute magnitude u'^ and points 
away from the Sun; thus i<"(^) has to point towards it. 
This is why LI is closer to the Sun than the SP of the 
potential. With the extra force (|81l) to balance, LI is 
further shifted toward the Sun. In view of Eq. JHUI) , the 
predicted shift is approximately 4 m. 

There is a similar order of magnitude effect on the or- 
bits about LI, and while this is not the primary purpose 
of the LPF mission, we suggest that a careful monitoring 
of the spacecraft trajectory may be of interest to gravi- 
tational physics. 



IX. THE VIEW INSIDE THE SPACE CAPSULE 

The form of the field (j) due to the Sun and Earth as 
discussed in previous sections is relevant for computing 
the orbit of a spacecraft or the effects of tidal stresses 
on experiments within it. Other questions are germane 
if we are interested in gravitational fields created by the 
spacecraft's components, as in a space reenactment of 
the Cavendish experiment. These are the subject of the 
present section. 

As long as the spacecraft's propulsion is off, it will 
move on a geodesic of the physical metric gafs of TeVeS 
basically because the energy momentum tensor of mat- 
ter is conserved with respect to that metric. If we ignore 
relativistic corrections, this path corresponds to a New- 
tonian trajectory in the potential <I>jv + 4' created by the 



SS. What we wish to ask is, in such an orbit what fields 
and forces exist within the spacecraft? 

If we proceed by analogy with general relativity, one 
might guess that in the spacecraft's frame of reference 
the sum of the perturbations S^j\j and (50 to <I>jv and 0, 
respectively, sourced by the spacecraft's structures and 
free proof masses would constitute the effective gravi- 
tational potential determining relative accelerations, etc. 
Below we show that this is so to sufficient approximation. 
In effect this result shows that TeVeS complies with the 
weak equivalence principle (in a freely falling frame ex- 
ternal gravity is cancelled out). TeVeS does not obey 
the strong equivalence principle. 

In TeVeS the physical metric (jap and Einstein metric 
ga/3 are related by (our signature is {— , +, +, +}) 



9ap 



21laH/3sinh(20), 



(82) 



where iXa is the eponymous vector field of the theory; it 
obeys ga/jil^il'^ = —1. Let gap and gap represent the 
metrics generated by the SS. In the frame of the freely 
moving spacecraft (supposed to be non-rotating) the 
physical metric induced by the SS will be of Minkowski 
form. The transition from the global frame to the space- 
craft frame is effected by projecting the said metrics with 
the help of a suitable tetrad; its explicit form will not be 
needed here. We use indeces a, 6, c, • • • to label spacecraft 
frame vectors and tensors. 

Now (j) of the SS is small; we regard its value within 
the spacecraft as a fixed number Likewise, Hq in 

the SS is a unit vector which points solely in the time 
direction (if we ignore the motion of the SS itself with 
respect to the Galaxy, etc.). We thus regard ita within 

the spacecraft as a fixed vector, ili"' , which in general has 
small space components of order the craft's velocity. We 
thus have for the Einstein metric within the spacecraft, 
to first order in (po and the spatial components of ili"-* , 



(83) 



The final term is the perturbation to the Einstein metric 
from the energy momentum tensor of spacecraft compo- 
nents. Within the spacecraft hab is the only part of gat 
whose space variation is significant. 

The gab metric comes from Einstein equations, as mod- 
ified in TeVeS ; here we work in linear approximation. 
As well known, the Einstein tensor can be linearized in 
terms of second derivatives of the perturbation to r/ab', 
only hab enters into it, in the form customary from gen- 
eral relativity, because the rest of the terms are here 
regarded as constant. Further, here as elsewhere, any 
raising of indeces can be done with 77°^ since the differ- 
ence between it and g"'' is already of first order in small 
quantities. The sources of the modified Einstein equa- 
tions contain scalar, vector and matter contributions. 
Most are quadratic in and ila derivatives, and thus 
quadratic in the small S(j) and 5iia corrections to cj)^^^ and 
Uq produced by the spacecraft. We can thus ignore these 
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energy-momentum contributions. Some further inspec- 
tion reveals that two other contributions, that related to 
the TeVeS Lagrange multiplier A, and that coming from 
the free function F, are likewise negligible. 

We conclude that the only source of hat is the mat- 
ter of the spacecraft. In TeVeS apart from the usual 
source Tab there is one of form Taf,(l — e'*'^). Obviously 
because of the smallness of this last is negligible. The 
linearized Einstein equations thus look like those in gen- 
eral relativity, and the relevant solutions for hab are the 
familiar ones. In particular, the temporal-temporal com- 
ponent is htt = —26^N- We compute the physical metric 
in the spacecraft's frame by substituting Eq. H83I) into 
Eq. ^ and setting e"^-^ = 1 - 2(0^") + (50) + • • • , 

sin(2(/)) = 2(0(°)) + (5( 
get 



• • • and ila = + SUa. We 



Vab + hab-2S(^ (?7„fc-f2ilWii[°^) 



(84) 



where the terms omitted are of second order in the small 
quantities S(j), Siia and hab- 

In calculating gtt we take note of the fact that 
differs from unity by a term of order the square of the 
spacecraft's velocity, which is of the same order as $Ar. 
Such a correction is negligible in Eq. We thus get 

gtt = — 1 — 2(5$jv + Accordingly, in the spacecraft's 
frame the physical gravitational potential equals 6^n + 
d(j) as surmised earlier. 

Now we know that S^n comes from Poisson's equa- 
tion. But how is S(p to be calculated? Let us substitute 
(f) = 0(0) -\- S(f> in the scalar equation Q . In the present 
calculations we shall take cognizance of the spatial gra- 
dient of 0^°-' and regard it as large compared to V(50. 
Linearizing the equation in derivatives of S4> leads to 



kGp, (85) 
(86) 
(87) 



AS(j) +2^W Wd^S(j)/dx' dx^ + ■ 
H = (|V0(")n-i/2v^(o)^ 

^ = dln^(y)/lny, 

where A represents the Laplacian and Y = fcP|V(/ 
In Eq. f H85|l the ellipsis denotes terms with first deriva- 
tives of 5(j) only. Quite in analogy with the eikonal ap- 
proximation we shall ignore these last; presumably the 
spacecraft's small scale makes contributions containing 
only a first derivative of a varying quantity subdominant. 

It should be evident that the unit vector H is antipar- 
allel to the U (discussed in Sees. IIIIVI|) coming from the 
SS. Aligning the coordinate system in the spacecraft's 
frame with its x axis in the H direction (possible at a 
particular position in the orbit), we see that Eq. H85|) 
is just a Poisson equation whose x coordinate has been 
rescaled to x(l + 2^)^^/^. Now for our model of 
n{Y) we calculate 



21+714 



(88) 



It follows that when the spacecraft is in the deep MOND 
region (/i built with V^'-'^-' is small compared to unity). 



we find the x direction is compressed by a factor 2. 
This is a facet of the "external field effect" whereby 
MOND effects in a weak field systems are traded for 
quasi-Newtonian behavior but with rescaling in one di- 
rection. 

By contrast, with the spacecraft deep in the quasi- 
Newtonian regime (^ « 1), 50 is determined by the 
usual Poisson equation. With (50 proportional to 5^m 
we only have a (small) rescaling of the effective gravita- 
tional constant, or equivalently of Milgrom's jl, a point 
which has already been mentioned in Sec. ^] Thus in a 
quasi- Newtonian environment, even the superweak fields 
originating in the spacecraft components behave in ev- 
eryday (Newtonian) fashion. 



X. CONCLUSIONS 

Postulating dark matter or MONDifying the gravita- 
tional interaction are conceptually conflicting ways of 
dealing with several anomalous astrophysical observa- 
tions. While it is possible that these anomalies will 
themselves decide between the two approaches, a "di- 
rect detection" would be far more convincing, for exam- 
ple, ongoing dark matter searches finding a particle with 
suitable cosmological and astrophysical features |23|. In 
this paper we examined what might constitute "direct" 
detection of MOND behavior. We predicted the exis- 
tence of regions displaying full MOND behavior well in- 
side the Solar system, specifically in bubbles surround- 
ing the saddle points of the gravitational potential. If 
abnormally high tidal stresses are observed in these re- 
gions this would prove MOND beyond reasonable doubt. 
Occasional astrophysical difficulties with the theory (e.g. 
with regards to lensing 21] or cosmological density fluc- 
tuations jl^ J^, ^] ) would no doubt mysteriously dissi- 
pate should such a discovery be made. 

How general are our predictions? MOND's solid re- 
quirement is that jjL{x) must approach 1 as x S> 1 and 
a; as a; ^ 1; the interpolating regime between these two 
asymptotic requirements is far less constrained. In the 
present work this intermediate regime translates into the 
quasi-Newtonian calculations presented in Sec. [V] For 
these we chose a reasonable form for ^jl{x), Eq. but 
we should stress that the details are model dependent. 
For instance, in H18I) the leading correction could have 
been quartic in uq/F instead of quadratic, resulting in 
a different power in the denominator of (|42(l . The extra 
force SF would then fall off more steeply with r. 

Accordingly, our calculations in the quasi-Newtonian 
domain are simply illustrative. We defer to a future pub- 
lication a thorough study of the effect of the choice of 
fi (as dictated by theoretical requirements and extant 
observations) on planetary orbits Lagrange points, 
and the Pioneer anomaly [9| . By contrast our predictions 
for the interior of the ellipsoid H38|) , as presented in Sec- 
tion IVII are robust predictions of the MOND scenario, 
and of wider validity. 
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We thus face a dilemma. The strongest MOND ef- 
fect and the theoreticaUy more robust prediction is that 
made in Section IvTl for the interior of the elhpsoid (j2HJ- 
However locating it in space may be taxing, particu- 
larly since this bubble is non-inertial. In contrast, the 
quasi-Newtonian predictions, e.g. what LISA Pathfinder 
might find in the vicinity of LI, are geographically less 
demanding, but the predicted effects are weaker and the- 
oretically less discriminative. Thus observing what we 
predicted in Section IVI would support the specific model 
(|24|l there; however failure to observe it would hardly 
disprove MOND in general. The interior of the ellipsoid 
(|38|l is therefore the prime experimental target for a con- 
clusive test. But one should not despair: systems other 
than those examined here may naturally reveal the inner 
core derived in Section IVII For example the movement 
of the saddle point through a diffuse medium - say the 
rings of Saturn - could be observable. 

On a technical note, our calculations once more un- 
derline the limitations of the usual folklore that MOND 
can be obtained by "taking the square root of the New- 
tonian field when V<I> <C ao". This is only true under 
strict spherical symmetry; in general one must add to 



the Newtonian field a curl term, which acts as a sort of 
"magnetic" gravitational field. As we have seen in the 
study of saddles, this field cannot be neglected even in 
the quasi-Newtonian region; and in the deep MOND re- 
gion neglect of the curl component introduces downright 
errors in the quantitative details. This should serve as a 
warning when making naive comparisons between theory 
and observation, for example in the context of clusters or 
satellite galaxies. 

There are other regions in the solar system where gra- 
dients of the Newtonian potential will be low, e.g at the 
center of near-spherical objects. However these are ob- 
viously inaccessible. By focusing on the saddle points of 
the gravitational potential in the solar system we believe 
we have exposed the best candidates for a direct detec- 
tion of strong MONDian behavior in our own backyard. 
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